5u(l,2) Algebraic Structure of the XYZ Antiferromagnetic Model in 

Linear Spin- Wave Frame 

Shuo Jin 1 *, Bing-Hao Xie 2 , Hong-Biao Zhang 3 , Jing-Min Hou 4 
1 Department of Physics, School of Science, 
Beihang University, Beijing, 100083, P. R. China 
2 Beijing Information Technology Institute, Beijing, 100101, P. R. China 
^ Institute of Theoretical Physics, Northeast Normal University, Changchun, 130024, P-R-China 
^Department of physics, Southeast University, Nanjing, 210096, P. R. China 

Abstract 

The XYZ antiferromagnetic model in linear spin-wave frame is shown explicitly to have an su(l, 2) 
algebraic structure: the Hamiltonian can be written as a linear function of the su(l,2) algebra gener- 
ators. Based on it, the energy eigenvalues are obtained by making use of the similar transformations, 
and the algebraic diagonalization method is investigated. Some numerical solutions are given, and the 
results indicate that only one group solution could be accepted in physics. 
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I. INTRODUCTION 



As is well known, the Heisenberg model is a simple but realistic and extensively studied 
solid-state system . According to the sign of the interaction intensity J, the model can be 
classified as the ferromagnetic type and the antiferromagnetic one. Based on the interaction 
intensity along the different space directions, the model can be labelled as the XXX, XX Z, 
and XYZ ones. Recently, it has been found that the Heisenberg interaction is not localized in 
spin system, and it can be realized in quantum dots [3, |4( , nuclear spins 0, cavity QED 
Thus, the study of this basic model is of interests and wide applications in physical fields. In 
the investigation of these models, it is a basic task to get the exact solutions 
i.e., to diagonalize the Hamiltonian. So far, only some special Heisenberg models can be 
exactly solved, such as XXX antiferromagnetic model ["L^ . In general, the linear spin- wave 
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161 ] approximation is widely applied when we study the Heisenberg models. It has been 
known that the XX Z model in linear spin-wave frame can be diagonalized by coherent state 

fin 

operators |17l.ll8| of su(l, 1) algebra. However, for the XYZ model in this frame, the method of 
the coherent states does not work. Therefore, it is necessary for us to develop another algebraic 
diagonalization method to get the energy spectrum by using the algebraic structure which the 
model has. 

In this Letter, we review the su(l, 1) coherent states of the XX Z antiferromagnetic model in 
linear spin-wave frame. Then, the XYZ antiferromagnetic model in linear spin-wave frame can 
be written as the generators of su(l, 2) algebra. At last, the energy eigenvalues are obtained 
in terms of the algebraic diagonalization method, and some numerical solutions are given and 
discussed. 

II. XXZ ANTIFERROMAGNETIC MODEL AND su(l, 1) COHERENT STATES 

The Hamiltonian of the XXZ antiferromagnetic model reads: 

Hxxz = - J E ( S ? S j + S f S ] + VS*S*) (J < 0), (1) 

<i,3> 

where the notation < i,j > denote the nearest neig hbor bonds. Starting from the two-sublattice 
model and Holstein-Primakoff transformation 

S z a = -s + a! a, S z b =s- 6 + 6, 
Si = (2 a )i(l-oto/2 S )ia, S- = (SlY, 



2 



Si = (2 S )i6t(i_ 6 t 6 / 2s )i i S^ = (Sl)\ (2) 

where a + and a, b) can be regarded as the creation and annihilation operators of boson 
on sublattice A (sublattice B), respectively, but the particle numbers n a = a + a, rib = b + b can't 
excel 2s, respectively. Because in low temperature and low excitation condition, < a) a >, < 
b^b s, the non-linear interaction in Eq. (J2J) could be reasonable ignored j^j]. Based on it, 
transferring the operators into momentum space, we obtain 

H xxz = -2ZsJ(Nsri-Y,H*), (3) 

k 

H k = r)(ala k + fc k 6 k ) + 7 k (a k 6 k + alMj- ( 4 ) 

Here 

7k = ^- 1 Ee ik ' R = 7-k, (5) 

R 

in which R is a vector connecting an atom with its nearest neighbor, and the sum runs over 
the Z nearest neighbors. 2N is the total number of the lattices. if k can be expressed as the 
linear combination of the su(l, 1) algebra generators in the form 

H k = 2 V E* + lll (E* + E k ), (6) 

with 

= E£ = aJ k , £ z k = -W + 4 + 1), (7) 

which obey the commutation relations of su(l, 1) Lie algebra: 

[E i l,E k ] = -2El [E*,E*] = ±El. (8) 

By introducing an su(l, 1) displacement operator 

W(&) = exp(&*$ - &E*) (9) 

with the coherent parameter £ k — r e* e , then we have 

W-\&)H k W(&) = aE* + ((3El + /T£ k ), (10) 

where 

a = 2r?cosh2r + 7k (e ie + e~ id ) sinh2r, (11) 
(3 = r)e ie sinh2r + 7 k (cosh 2 r + e 2ie sinh 2 r). (12) 
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In order to diagonalize H^, the coefficient before the non-Cartan generators and of Lie 
algebra, /3, should be chosen to zero (here we set 9 = for simplicity) and this leads to 

tanh2r = -^, a = 2yV - 7k \ (13) 

So if we denote 

\&>=W{&)\vac), (14) 

then one has 

#k|£k >= (n a + n b + l)e k |f k >, (15) 

e k = 2ZJsy/rP->£, (16) 
where requiring |7y| > |7 k |. The diagonalization of the XX Z antiferromagnetic model in linear 
spin-wave frame has turned out to be a direct product of su(l, 1) coherent states (8> k |£k)- One 
can see e k is the quantum of antiferromagnetic spin-wave, i.e. the dispersion relation. From 
Eq. ()15|) . it is known that for any k, there exist two branches of degenerate antiferromagnetic 
spin-wave which the quasi-particle numbers are described by n a and n b , respectively. 



III. XYZ ANTIFERROMAGNETIC MODEL IN LINEAR SPIN-WAVE FRAME 
WITH THE au(l,2) ALGEBRAIC STRUCTURE 

Owing to the different interaction intensity along the different space directions, in general, 
the Hamiltonian of the XYZ antiferromagnetic model is described by 

H XYZ = -J E {v x S?Sf + r] y S?S? + S?S*) 

<i,j> 

(J<0, 7fc, % >0), (17) 

where we have set r\ z — 1. Similar to the former case of the XX Z antiferromagnetic model, 
the XYZ antiferromagnetic model in linear spin-wave frame is given by the Hamiltonian: 

H XYZ = 2ZsJ[Ns - (]Tft k - 1)], (18) 

k 

H k = a k a k + 6 k &k 

+t; k (a k &L k + 4.6-k) + Pk(a k &k + albl), 

(19) 



with 



v, = ^7,, = (20) 
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If we choose 



/+ = a£6_ k , / k = a k &l k , = a k 6 k , 
V* = b+b\, V k = M>_ k , U* = a+b+, 

/ 8 k = 4K + ^-k + 24 + 2), 



(21) 



then, one can see they obey the commutation relations of sw(l,2) Lie algebra (here we omit 
the momentum sign k): 

[J 3 , J±] = ±J±, [/+, J_] = 2/ 3 , [Is, 4] = 0, (a = ±, 3) 
[/ 3 , C4] = T^±, [4, E4] = ±U ± , [U+, C/_] = I 3 - ^s, 
[/ 3 , 14] = t\v±, [Is, V ± ] = TV ± , [V + , V.} = 7 3 + ~/ 8 , 



[J ± , C/ ± ] = T V T , [U±, V ± ] = ±I T , [7±, F ± ] = ±f/ T , 



[J±,l/ T ] = [J ±) V T ] = [17 ± ,V T 



0. 



(22) 



From Eqs. (fK?j) and (j21j) . 7i k can be expressed as the linear combination of six generators of 
Lie algebra su (1,2) and processes su(l, 2) algebraic structure, i.e., 



rk 3 rk 



/ 8 k +Pk(/:+/ k )+^(f/ k + f/ k : 



(23) 



IV. THE DIAGNOLIZATION AND THE EIGENVALUES 



If we set the general linear combination form of Lie algebra su(l, 2) as 



H = al + + &/_ + cU+ + d£7_ 



(24) 



for 7i k the coefficients in Eq. (|2*4*jl are: 



a = b = p k , c = d = t?k, e = / = 0, 
3 

9 = 1, h = ~^- 



(25) 



Until now, the coherent state operator as Eq. Q has not been found for the XYZ antifer- 
romagnetic model in linear spin-wave frame. But following the standard Lie algebraic theory 
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24j, if H is a linear function of the generators of a compact semi-simple Lie 



group, it can be transformed into a linear combination of the Cartan operators of the corre- 
sponding Lie algebra by 

Hi = WH W~. (26) 

Here W = lli=i ^xp(xiAi) is an element of the group and W~ denotes the inverse of W, in 
which Ai (i = 1, ...,N) is a basis set in Cartan standard form of the semi-simple Lie algebra, 
and Xi can be set to zero if the corresponding Ai is a Cartan operator. By choosing 

W = exp(x 3 iV + )exp(x 2 ilJ)exp(x 32 UJ) 

exp(x 12 I + )exp(x 23 U + )exp(x 13 V^), (27) 

and letting the coefficients of the non-Cartan operators vanish, while substituting Eqs. 
(j2U) (jlil) (J2Zj) into the right-hand side of Eq. (}2T)j) . we get a complete set of algebraic equa- 
tions of Xij after lengthy computation: 

-(h + lg)x 13 + (a + dx l3 )x 23 = 
' c + bx 13 + (\g - h)x 23 + dx 2 23 = (28) 
a + dxi 3 — (g + dx 23 )x\2 — bx\ 2 = 0, 

and the Hamiltonian after the transformation of W becomes diagonal: 

Hi = WH W- 

3 

= (g + dx 23 + 2bx l2 )I 3 + (h - -dx 23 )I 8 . (29) 

One can see that although operator W is not unitary, the similar transformation (}2l)j) guarantees 
that the eigenvalues of Ho equal those of H±. This is acceptable for we are only concerned 
with the eigenvalues. If the total particle number Af = n£ + n b _ k + n£. (n^ = a^a^, n'L k = 
fr-k^-k, K = fej^k), from Eq. (J21, [T,M] = (T = J±, V±, U±, I 3 , I s ) holds. Hence, supposing 
the common eigenstates of the Cartan generators -Z3 and of Lie algebra su(l, 2) are the Fock 
states I n^_,n b _ k ,n^ >, i.e., for the commutative set {I 3 ,I$,Af} there exist: 

h I <,n\,nl >= \{K-n\) \ <,n b _ k ,4 >, 
h I K,n b _ k ,n b k > 

= -~« + n\ + 24 + 2) I n a k , n b _ k , 4 >, 
M I K,n b _ k ,ni > 

= K + n b _ k + n b k ) \nln b _ k ,n k > . (30) 
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From Eqs. (|29|) (|30 |) it follows the eigenvalue of the Hamiltonian (J23j) : 



UJ 




(31) 



where the coefficients b,d,g,h, are given in Eq. (J25j) and x^ can be obtained by solving Eq. 
(f2"%J) . and u;^, u; k , k are the energies of the three different magnons respectively. In fact, 
the order of the operators in W can be chosen arbitrarily, but the coefficients Xi are strongly 
dependent on the order. Although any specified order has a solution, a properly chosen order 
can simplify the procedure to get the X{. In general, for the Hamiltonian with su(n) (whose 
Cartan operators are An = bfbi) or the isomorphic algebra su(p, q) (p + q = n) structure, the 
transformation operator W can be chosen as 



where the order of the operators of exp(xijAij) (i ^ j) is arranged according to the roots of A^- 
in a decreasing way. For example, the root of A^i is highest, and that of A^ is lowest. As a 
rule of our choice, the right or left Aij in Eq. ()32j) is the one that is missing in the Hamiltonian 
7io; the middle operators sequence forms a circle root diagram. With this specification, in our 
experience the coefficients Xij are relatively easy to work out. We also choose the form of W in 
Eq. (J27)) as exp(xi3V-)exp(xi2l+)exp(x23U + )exp(x2iI-)exp(x32U-)exp(x3iV + ), and it can be 
proofed they lead to the same eigenvalues. 

V. NUMERICAL SOLUTIONS 

From the Eq. (|2"8j) . maybe there have several sets of solutions, which consist the complete 
set. But the different sets of the slutions cannot be the eigenstate of the Hamiltonian (|19|1 or 
()23j) simultaneously. Only those who possess the physical meaning is the solutions we need. 
In order to illustrate it, we consider a concrete example of the Hamiltonian (|19|) or (J23|) with 
r] x = 0.8, r/y = 0.5,7k = 1. Then t> k = 0.15, p k = 0.65. Using maple, one can show that there 



W = exp(x N iA N1 )exp(x N2 A N 2)... 



exp(x 2N A2N)exp(x 1N A 1N ), 



(32) 
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are six sets solutions of Eq. ([28)1. in which only one set x i3 = — 0.6018692595e — l,x 2 3 = 
0.9389946768e - l,x 12 = 0.4827143955 leads to the positive energy = 1.327849277, u;£ = 
1.014084920, u\ = 0.6862356430 of the magnons. It is clear that only this solution is the 
accepted in physics. Other five sets solutions (with negative energy) are non-physical. This 
procedure is easily to be taken in solving the Eq. (|2"%j) and Eq. ((221) • 

VI. CONCLUSION AND REMARKS 

In conclusion, the eigenvalue problems of the XYZ antiferromagnetic model in linear spin- 
wave frame is solved from an algebraic point of view. To use this algebraic diagonalization 
method, first, it is needed to find the algebraic structure of the Hamiltonian and manage to 
write the Hamiltonian into a linear combination form of the algebraic generators just as Eq. 
()23|) . Second, according the particular structure of some Lie algebra, one may looking for 
the transformation operator. The key is that we can let the coefficients of the non-Cartan 
operators vanish successfully and get the solvable equations of the parameters through the 
transformation. Some numerical solutions check our diagonalization method, whose advantage 
is that the eigenvalues of Ho equal those of Hi although the Hamiltonian changes. Due to 
the different interaction intensity along the different space directions for the Heisenberg model, 
the complication in physical model leads to the enlargement of the algebra structure, such as 
the XX Z case to the XYZ case corresponds to the su(l, 1) algebra to the sm(1,2) one. Of 
course, the change of the algebra structure brings the different method in diagonalizing the 
Hamiltonian. It is reasonable to believe that more useful physical applications of the algebraic 
diagonalization method should be found. It may be possible to extend this case to higher-rank 
Lie algebras. 
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